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We analyze the number of particle fluctuations within the semiclassical approximation of an ideal 
bosonic gas with an anomalous single-particle dispersion relation suggested in several quantum- 
gravity approaches trapped in a generic 3-dimensional power-law potential as a criterium of ther- 
modynamic stability for these systems. We show that the analysis of the fluctuations in the number 
of particles caused by a deformation in the dispersion relation, leads to an observable consequence 
associated to the thermodynamic stability for a specific choice of the trap parameters. Additionally, 
"^j". we deduce the shift in the critical temperature in the thermodynamic limit associated with this 

"modified bosonic gas" and show that this shift expressed in function of the number of particles, 
can be used as an amplifier for some quantum gravity manifestations for stable systems. 
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I. INTRODUCTION 



The search for small manifestations of quantum gravity or "quantum gravity phenomenology" in our low energy 
world is a very controversial topic in modern physics. In some schemes, the possibility that the space-time could be 
quantized, can be characterized, from a phenomenological point of view as a modification in the dispersion relation of 
microscopic particles [l|-y| • A modified dispersion relation emerges as an adequate tool in the search for phenomeno- 
logical consequences caused by this type of quantum gravity models. Nevertheless, the principal problem in the 
search of manifestations of quantum gravity is the smallness in the predicted effects [3j . If this kind of deformations 
are characterized by some Planck scale, then the quantum gravity effects become very small [HQ, but in principle, 
different from zero. In the non-relativistic limit, the dispersion relation can be expressed as follows Q 



E ~ m 
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Equation (Q} is expressed in units with the speed of light c = 1, and M p (~ 1.2 x 10 28 ey) is the Planck mass. The 
three parameters £1, £2, and £3, depend of the quantum gravity model in question [3, [H, and should take positive 
■ or negative values close to 1. Equation ([IJ, is the starting point in the search for small manifestations of quantum 
I gravity at low energies. 

The Bose-Einstein condensation phenomenon, from the theoretical and experimental point of view has produced 
an enormous amount of publications associated to this topic (3-l2(jj. Amon g th e issues addressed we may find, its 
. possible use as tools in the search of quantum-gravity manifestations [H, l2lM25| . For this propose, let us define the 
next "modified Hamiltonian" 



P 2 

H = i— + ap + U{r). (2) 
2m 

Where p is the momentum, m is the mass of the particle, and the term ap, with a = £1 ^g-c in ordinary units, is 
the leading order modification in expression ([1]), being c the speed of light. 
The potential term 



U(r) = J2A^ , (3) 
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in the Hamiltonian ([2]) is the so-called generic 3-dimensional power-law potential, where Aj and a, are energy and 
length scales associated to the trap 13]. On the other hand, fj are the d radial coordinates in the n^-dimensional 
subspace of the 3-dimensional space. The sub-dimensions rij satisfy the following expression in three spatial dimen- 



sions 



(4) 



If in equation Q d = 3, n\ = 7J2 = n% = 1, then the potential becomes in the so-called Cartesian trap. If d = 2, 
m = 2 and ri2 = 1, then we obtain the cylindrical trap. If d = 1, ni = 3, then we have the spherical trap. On the 
other hand, if s; — > oo, we have a free gas in a box. In this sense the potential included in the Hamiltonian @ is 
quite general. Different combinations of these parameters give different classes of potentials, according to ([3]). 

The study of fluctuations in the number of particles in the Bose-Einstein condensation phenomenon is a very 
important topic [l5l-l20| . First of all, the fluctuations in the number of particles are directly related with the equivalence 
or non-equivalence of the statistical ensembles which is a non-trivial and deep topic [26| . On the other hand, the 
fluctuations in the number of particles are directly related with the thermodynamic stability of the system in question 
[lil [l9l | . The fluctuations in a Bose-Einstein condensate depend strongly on the single-particle energy spectrum 
of the trap [20| and, consequently, the modification in the dispersion relation (fTJ) shall affect the thermodynamic 
stability of the condensate, and could lead to an observable manifestation with an appropriate manipulation of the 
trap parameters, defined in the generic potential ([3]). 

In this work, we analyze the Bose-Einstein condensation phenomenon within the formalism of the semiclassical 
approximation, assuming as a fundamental fact a deformation in the dispersion relation. Afterwards, we analyze 
the effects of this assumption upon the thermodynamics of the considered bosonic gas, explicitly, over the critical 
temperature. Additionally, we study the thermodynamic stability of these systems through the fluctuations caused 
by the deformed dispersion relation (JTJ) . Finally, we give the main results and adds some comments about the possible 
detection of some manifestations of quantum-gravity using different kinds of potentials. 

II. SEMICLASSICAL APPROXIMATION FOR AN IDEAL MODIFIED BOSONIC GAS 

The first part of the present work addresses the issue of the effects of a deformed dispersion relation over the 
condensation temperature for an ideal bosonic gas, trapped by the generic potential given by expression ([3]) in the 
semiclassical approximation 0, H3| • 

The semiclassical energy associated with the modified Hamiltonian @ is given by 

2 

e i^P)=n ^ap+U(r). (5) 

In the semiclassical approximation, the single-particle phase-space distribution may be written as [4, 27] 

n ^ = e mr, P i-^-i ^ 

Where (3 = 1/kT, being k the Boltzmann constant and T the temperature. Additionally,/^ is the chemical potential. 
The number of particles in the 3-dimensional space obey the normalization condition [J, [2^] , 



where 



N = I d s rd s pn{f^ = I d s fn(r) = / d'pnip), (7) 



n(r) = J d 3 pn{f,p), (8) 
n(p) = [ dVn(f, P | (9) 
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are the spatial and momentum densities, respectively [H1I3]- Using the expression ((5|), and integrating expression 
© over the momentum space, using expression©, allows us to obtain the spatial distribution associated with the 
modified semi-classical spectrum §5§ 



n(r) 



where 



(10) 
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is the dc Broglie thermal wavelength, jU e // = A 1 + ma 2 /2, and g v {z) is the Bose-Einstein function defined by [26 



9u{z) 



x v ~ x dx 



Y(v) J z~ x e x -\' 



(12) 



Where z is the so-called fugacity, which is related to the chemical potential through z = e' 3 ^ [26], [27j ■ The Bose- 
Einstein function (fT2j) diverges for z = 1 when i/ < 1 [26| . 

If we set a = in equation (|10j) we recover the usual result for the spatial density in the semiclassical approximation 

0. EE}. 



(13) 



Integrating the normalization condition ([7]), using expression (ITC)|) . with the corresponding potential ([3]), allows us 
to obtain an expression for the number of particles in function of the chemical potential fi, the temperature T, and 
the deformation parameter a, given by 



N 



N a + CHf =1 A, a^TQ + l) (-^) 3/2 57 ( Ze//)(KlT 
a \2^) 9 ^~ 1/2 ^ eff ^ KT ^ 1 2 \2^) f 7 -i( z e//)(«:T) 7 - 1 



(14) 



Where 



~ 2 ^ si ' 



(15) 



is the parameter that defines the shape of the potential. Nq are the particles in the ground state, T(y) is the Gamma 
function, and C is a constant associated to the potential in question. In the case of Cartesian traps, and in consequence, 
for a three dimensional harmonic oscillator potential C = 8. Additionally, z e ff = e Pd J -+ ma / 2 ) = ze ma / 2kT is an 
effective fugacity. Setting a = we recover the usual expression for the fugacity z [26j. 

Accepting that ma 2 /2 << kT, we can expand the expression (fT4l around a — to first order, using the properties 
of the Bose-Einstein functions I2a| 



(16) 



Using this fact we can re-write (fT4"f as follows 
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Let us define 



(17) 



(18) 



as the characteristic volume associated with the system . We can see that from expression ([18)1 . when Sj — >• oo, then 
V c / la r becomes the volume associated with a free gas in a box. In this sense Vchar can be interpreted as the available 
volume occupied by the gas [IS Ell- O n the other hand, it is noteworthy to mention that the most general definition 
of thermodynamic limit can be expressed as 



N -> OO, V cha r -> 0, (19) 

keeping the product NV c har — * const, and is valid for all power law potentials in any spatial dimensionality (l9j . 
With the criterium given above, the critical temperature in the thermodynamic limit is well defined. 

At the critical temperature in the thermodynamic limit fi = and No = 0, which implies that the Bose-Einstein 
functions are given by the corresponding Riemann Zeta functions £(x) [26| . then from expression (|17p and the definition 
([T5| we obtain that the number of particles at the condensation temperature can be written as 



N 



V c ha 
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(20) 



where T c is the critical temperature for a modified bosonic gas in the thermodynamic limit. It is noteworthy to 
mention that the phenomenon of condensation for a modified bosonic gas in the thermodynamic limit is possible when 
7 > 3/2 according to the properties of the Bose-Einstein functions (IT2)) 26]. When s, — > oo, 7 = 3/2, we have a free 
gas in a box, in this case apparently the Bose-Einstein condensation is not possible for a free modified gas in a box, 
because the divergent behavior of the Bose-Einstein functions ([12")) . To make possible the condensation of a modified 
free gas in a box we have to take into account the minimal energy associated to this system [25[ . 

From expression ()20[) . setting a — 0, we may obtain the expression for the critical temperature To for a gas trapped 
in a generic 3-dim power-law potential in the thermodynamic limit [l3| 



To 



N 
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From equations ([20)l and (f2~T)) we obtain the shift in the critical temperature for our deformed bosonic gas 



(21) 



H 7 =M 7+a (v) 



2m\ !/2£( 7 - 1/2) 



C(7) 



kT 



7-1/2 



(22) 



We can notice from expression (|22[) , that T c increases for a > 0, compared with the corresponding critical tem- 
perature of the usual Bosonic gas when a = 0. The opposite case, a < 0, corresponds to a decrease in the critical 
temperature, compared with the usual case. The effect of the external potential is to concentrate particles in the 
center of the trap and a positive a increases the critical temperature which means that the effect of a > "reinforce" 
the main effect of the external potential. When a < we have the opposite behavior, a negative a tends to "weaken" 
the effect of the external potential. We can also express the shift in the critical temperature as a function of the 
number of particles N, 



To 
Where 



A/ * a nN^/ 2 \ (23) 
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= (2m\ 1/2 C( 7 -l/2) / (2mh^ \ ^ 
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It is noteworthy to mention that the correction in the critical temperature depends strongly on the functional form 
between the number of particles and the parameters of the potential. This fact could be used as an amplifier of some 
quantum gravity manifestations in systems with finite number of particles [25j . 

The appearance of a in expression (1221) or (|23p may be interpreted as a parameter that modified the thermodynamic 
stability of the system. One way of quantify the thermodynamic stability of the system is through the fluctuations in 
the number of particles , which are directly related to the so-called isothermal compressibility [18|, [l9|, [26[ . 

Let us suppose now that our condensate is trapped in an anisotropic three-dimensional harmonic-oscillator potential. 
For this trap the the shape parameter is given by 7 = 3 with A% = hui/2, a, = y/h/rnuii (see expression (|18p). and 
using the definition <D = (wiCl^c^) 1 / 3 expression (|23p becomes 

ATc ~- Wr«i^\ (25) 



Ti 



~ 3C(3)5/^ 




III. FLUCTUATIONS IN A MODIFIED BOSONIC GAS 

The number of particle fluctuations are characterized by the dispersion [lH [Til [26| 



(AN) 2 = (^V 2 ) - (AO 2 = «r( (26) 

Where A" is the total number of particles of the gas. The over-line in expression (|26[) means an average over the 
ensemble, which corresponds to the thermodynamic variables [26]. When (A7V) 2 ~ A" the fluctuations are called 
normal, while if (AA^) 2 ~ N s the fluctuations are anomalous, with s > 1 [l8|, The analysis of fluctuations in 
statistical systems is important because the fluctuations define for instance, the stability of the system and how this 
system leads to the state of thermodynamic equilibrium [IH, . 

The expression (|2l)l) is directly related to the isothermal compressibility kt through the next relation [l8|, l26j 



(AN) 2 
pNnT ' 



(27) 



Where p is a mean particle density. Using expression (|I7[) . we may calculated the fluctuations associated with our 
modified bosonic gas in the thermodynamic limit. A necessary condition for the stability of a system is the semi- 
positiveness and finiteness of the isothermal compressibility. If the isothermal compressibility is negative or divergent, 
could lead to a immediate collapse or explosion of the system. In this aim, the analysis of fluctuations caused by the 
quantum structure of space-time in the number of particles, must lead to a observable consequences for a specific 
choice of the trap parameters. It is means, that the fluctuations caused by the deformation in the dispersion relation 
allows us, in principle, discriminate in which kind of systems the manifestations of quantum gravity are more feasible 
to be detected by analyzing the stability of the system. 

Assuming that N — N e , being N e the number of particles in the excited states, allows us to calculate from (|17[) . 
together with (|2l)|) , the fluctuations associated with our modofied bosonic gas 



(AAy 2 = [kT)W cHc 



wo /on 27r(2m) 3 / 2 Sirm 2 ^_i /2 , , 

F(3/2) (2^)3 9 ^ Z) ~ a (2^ (KT) 9 ^ {Z) 



(28) 



This last expression is a direct consequence of at = a^o + Ar e m m. From the relation between the critical 

temperature and the characteristic volume (|2I|) . allows us to write (|28|) above the critical temperature (T > To) as 
follows 
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(29) 



For T > Tq the fugacity < z < 1 [26[, and the Bose-Einstein functions in expression (1291) have well behavior. 
The appearance of a increase the fluctuations if a < 0. On the other hand, a positive a decreases the fluctuations. 
Nevertheless, the corrections caused by the deformation leaves the fluctuations normal, independent of the sign of 
a. The isothermal compressibility is always finite and positive for 7 > 1 and, hence, the system is always stable for 
temperatures above To. 

For temperatures T < T the dispersion ([29| becomes 



(a Ay 



2 _ 



C(7 - 1) 
C(7) 




1/2-1 




(30) 



From expression (|30[) . we notice that there is a critical value at which the fluctuations go from being normal to 
be anomalous, this value is 7 = 5/2. For values 7 < 5/2 the fluctuations becomes anomalous, and the isothermal 
compressibility becomes divergent and negative when a > 0, which means that the system is unstable. 

Setting a = in expressions (f2l?|) and (|3T)|) we recover the results given in (l8| . 

Let us analyze the case of an anisotropic harmonic oscillator in 3-dim. In this case 7 = 3, then from expression 
(f2"9"| . we obtain 



(AiVe 



52(2) 33/2(2) / 8m 



1/2- 



C(3) 



C(3) 
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(31) 



Expression ([31]) shows that the fluctuations are normal, independent of the sign of a. Then, the system is stable 
for temperatures T > Tq. On the other hand, for temperatures T < Tq, we obtain from expression ([30| 



(AN e y 



C(2) C(3/2) / 8m 



1/2- 



C(3) C(3) W«T 



N. 



(32) 



The fluctuations in the number of particles for temperatures T < Tq are normal, which means that the isothermal 
compressibility is finite and positive. Hence the system is stable beyond the critical temperature. We notice from 
expressions (|3"Tj) and ([3"2l that the modifications caused by a increases the fluctuations when a < 0, and decreases 
them in the opposite case i.e., a < 0. Nevertheless, the fluctuations are normal in both cases. The system is stable 
for temperatures above and beyond the critical temperature in the case of a harmonic oscillator potential. 

Setting a — in expression (I32p for temperatures T < T we recover the result given in (l5j . 



(AA^) 2 




02) 
C(3) 



N. 



(33) 



IV. CONCLUSIONS 



Using the formalism of the semiclassical approximation, we have analyzed the Bose-Einstein condensation for 
modified bosonic gas trapped in a 3-D power law potential. We notice that the critical temperature must be corrected 
as a consequence of the deformation in the dispersion relation. The shift in the critical temperature depends strongly 
on the sign of a. For a positive a, the critical temperature increases respect to the usual value when a = 0, and 
decreases in the opposite case, a < 0. The main effect of the trapped potential is to concentrate particles in the center 
of the trap. The insertion of the deformation parameter a modifies this behavior, depending on its sign. This fact 
allows us to interpret a as a parameter related to the thermodynamic stability of the system, since its appearance 
modify the fluctuations in the number of particles, and the isothermal compressibility, which is in an intimate relation 
with such stability. We notice that the fluctuations for T > Tq, remains normal for values 7 > 1 of the shape 
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parameter. On the other hand, a detailed analysis of the different values for the shape parameters, below the critical 
temperature T < T , shows that the system is stable for 7 > 5/2. A very important consequence from expressions 
(fTT|) and (|32|) is that for values between 3/2<7<5/2we have condensation, but the system is unstable. In this 
range of values of the shape parameter 7 we have, for example, potentials of the type V(r) ~ x 3 + y 3 + z 3 in the case 
of a Cartesian potential, V(r) ~ p 3 + z 3 in the case of cylindrical traps, and V(r) ~ r 3 in the case of spherical traps. 
We have noticed also, that for potentials of the type V(r) ~ p 4 + z 2 , V(r) ~ r 6 or any other combination such that 
3/2<7<5/2we have condensation, but the system is unstable. For systems with 7 < 3/2 (in which is included the 
ideal modified gas in a box) the condensation apparently is not possible if ones assumes that the minimal energy of 
the system is zero, or equivalently that the value of the chemical potential fi is zero at the transition temperature. 
Nevertheless, an inspection of the deformed relation dispersion (TTJ) show that the minimal energy is not zero. Then, 
in principle, the condensation is possible even in the case of a modified gas trapped in a box (2f| . In the case of stable 
systems 7 > 5/2, like the harmonic oscillator in 3-D, for example, one possible manner for the analysis of quantum 
gravity manifestations in a Bose-Einstein condensate is through the study of the thermodynamic variables like the 
critical temperature in function of the number of particles, as can be seen from expression ([22]) . Systems with a large 
but finite number of particles, in the ideal and weakly interacting case, can be used, in principle, as an amplifier of 
some quantum-gravity manifestations. A detailed analysis of these topics is analyzed in extenso in [25| . 
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